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Abstract. In this paper we describe an efficient involutive algoritlim for construct- 
ing Grobner bases of polynomial ideals. The algorithm is based on the concept of 
involutive monomial division which restricts the conventional division in a certain 
way. In the presented algorithm a reduced Grobner basis is the internally fixed sub- 
set of an involutive basis, and having computed the later, the former can be output 
without any extra computational costs. We also discuss some accounts of experi- 
mental superiority of the involutive algorithm over Buchberger's algorithm. 



1. Introduction 

The concept of Grobner bases for polynomial ideals together with the algorithmic char- 
acterization of these bases in terms of 5— polynomials was invented almost 40 years ago 
by Buchbergerin his PhD thesis 1 1 1. Since that time Grobner bases have become the most 
universal algorithmic tool in commutative algebra and algebraic geometry (see, for exam- 
ple, books |2|-|8 1). Moreover, this role of Grobner bases was extended to rings of linear 
differential and shift operators and to some other noncommutative algebras (see L9.10J 
and references therein). 

The fundamental property of a Grobner basis of a polynomial ideal is the divisi- 
bility of any element in the related initial ideal by the leading monomial of an element 
in the basis. In our paper II II an algorithmic approach was developed based on an ax- 
iomatically formulated restriction for the conventional monomial division. This restricted 
division which we called involutive generalizes the pioneering results by Zharkov and 
Bhnkov II12I on carrying over the involutive methods from differential equations to com- 
mutative algebra and leads to a more general concept of involutive bases. An involutive 
basis is a (generally redundant) Grobner basis. But for all that, an element in the initial 
ideal is multiple of the leading term of an element in the involutive basis not only in 
the conventional sense but also in the involutive sense. Just as a monic reduced Grobner 
basis 1131 . a monic minimal involutive basis L14J of an ideal is uniquely defined for a 
given monomial order. 

An impetus was given to the involutive division technique with the experimental 
demonstration in I I 11121 the efficiency of involutive algorithms as an alternative method 
to compute Grobner bases. The concepts, ideas and algorithms developed in LI IJ got 
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further development and extension to linear differential algebra and noncommutative 
rings in our papers fT4l - |[T7l and in the papers of other authors |[T8l - l26l . 
With all this going on, the involutive algorithms have been implemented for particular 
involutive divisions in Reduce 1,11.12 ,16,1 . Mathematica [2U, Maple MuPAD L29J . 
C/C++ lIT^ and recendy in Singular (version 2.0.5). 

Having analyzed our paper 1 1 1 1, Apel 1 30] introduced a slightly different concept of 
involutive division and designed another form of involutive algorithms based on his con- 
cept. The approach of Apel was further elaborated by himself and Hemmecke 1311 - 1341 . 
Differences of both approaches ITTl and BOl and some unifying aspects were studied 
in I.35J . 

For a given finite polynomial set, an involutive division partitions the variables into 
two disjoint subsets called multiplicative and non-multiplicative. The construction of an 
involutive basis is often called completion to involution and consists in examining the 
involutive reducibility of non-multiplicative prolongations, i.e., products of the polyno- 
mials and their non-multiplicative variables. In doing so, the involutive algorithms also 
examine all essential S*— polynomials |14 30| but with rather special reductions based 
on involutive division. Thereby, there is no conceptual distinction between Buchberger 
algorithm and involutive algorithms. The difference is in both critical pairs selection and 
their reduction provided with the presence of extra polynomials at involutive completion 
which are redundant in the conventional sense. 

Computer experiments with our C/C++ code 1 16 1 implementing the involutive algo- 
rithm for Janet division and with the standard benchmarks 1 36 37 1 show that in the most 
cases our code is faster than the most optimized modern-day implementations of Buch- 
berger's algorithm ' . Thus the involutive algorithm may be considered as an improve- 
ment of Buchberger's algorithm. 

In the present paper we follow the concept of involutive division in paper 1111 . which 
is briefly described in Sect. 3 as a particular case of the restricted monomial division 
introduced in Sect.2, and present our involutive algorithm in its detailed and optimized 
form (Sect.4). The algorithm can also be used to output the reduced Grobner bases of 
polynomial ideals without extra reduction of the computed involutive basis. Some aspects 
of computational efficiency of the involutive algorithm are considered in Sect. 5. It is easy 
to rewrite the algorithm for linear differential polynomial ideals and for ideals in rings 
of noncommutative polynomials with algebra of solvable type. Besides, the algorithm 
admits a straightforward extension to modules in these rings. 

2. Preliminaries 

In this paper we shall use the following notations and conventions: 

• N>o is the set of nonnegative integers. 

• X = {xi, . . . , is the set of polynomial variables. 

• X C X is a subset (possibly empty) of the variables. 

• K = K[X] is a polynomial ring over a zero characteristic field K. 

• Id(F) is the ideal in M generated by C K. 
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• M = {x^-l • • • a;^" \ik ^ N>o, 1 < A: < n} is the monoid of monomials in M. 
Mjc C M is a submonoid of M generated by the power products of variables in 
X. 

• degj {u) is the degree of Xi'mu E M. 

• deg(M) = dcgi(u) is the total degree of u. 

• is an admissible monomial order such that xi )^ X2 )^ ■ ■ ■ )^ Xn- 

• u \ V \s the conventional divisibiUty relation of monomial v by monomial u. If 
u I V and deg(u) < deg(z)), i.e. it is a proper divisor of v, we shall write u\Zv. 

• lm(/) and lt(/) are the leading monomial and the leading term of / G M \ {0}. 

• lm(F) is the leading monomial set for F c K \ {0}. 

• £ is an involutive monomial division. 

• u\c vis the divisibility relation of monomial v by monomial u defined by invo- 
lutive division C. 

• C{U) is the cone generated by monomial set /7 c M. 

• Cc{U) is the £— cone generated by monomial set J7 C M. 

• NMc{u, U) is the set of £— non-multiplicative variables of monomial u & U C 
M. 

• HNFc if, F) is the £-head normal form of / e M modulo F cR. 

• NFc (/, F) is the £-(full) normal form of / e M modulo F cR. 

• pol(p) = / is the first element in triple p = {/, g, X} where f,g€:FcR. 

• anc(p = {/, g, X}) = p is an ancestor of f in F. 

• nmp(p = {f,g,X}) = X is a subset of non-multiplicative variables for /. 

• idx(/, F) is the index of an element / G F in an indexed set F. 

To introduce the concept of involutive division and describe the involutive algorithm for 
construction of Grobner bases we need the following definitions. 

Definition 2.1. A linear monomial order >- is called admissible if the conditions 

m ^ 1 <^=^ m >- 1, nil >- ^2 <^=^ niim >- m2m 
hold for any monomials m, mi , m2 G M . 

Hereafter we consider only admissible monomial orders and omit the word "admissible". 

Definition 2.2. For a finite set C K \ {0}, a polynomial /i G K \ {0} and a monomial 
order >- consider polynomial h given by 

h = h -'^aijiTiijfj , 

ij 

where aij G K, fj G F, rriij G M and \m(mijfj) ^ \m{h) whenever aiij ^ 0. If 
lm(/i) has no divisors in Imi^, then h is called a head normal form of h modulo F. In 
this case we shall write h = HNF{h, F). In addition, if all other monomials occurring 
in h also have no divisors in \m{F), then h is called a normal form of h modulo F. In 
such a case we shall write h = NF{h, F). 

Definition 2.3. Given an ideal I c M and an order y, a finite subset G C M is called 
Grobner basis of X if 
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(V/eZ)(3geG) [lm(5) |lm(/)], 



(1) 



where u \ v denotes the conventional divisibihty of monomial v by monomial u. 
From Definition 2.2 it follows that Q is equivalent to 



Definitions 2.1-2.3 are well-known ingredients of the Grobner bases theory Q-O. Now 
instead of the conventional monomial division in ([0 we consider another divisibility 
relation and the related concepts of the normal form and Grobner basis. 

Definition 2.4. A restricted division r on M is a transitive relation u \r v (u,v e M) 
such that 



If u \r V, then u is r~divisor of v and v is r-multiple of u, respectively. Note, that the 
whole class of restricted divisions includes the conventional division as well. 

Definition 2.5. Let ;^ be a monomial order, r be a restricted division and F C M \ {0} 
be a finite polynomial set. If g e R \ {0} has a monomial which is r— multiple of an 
element in Im(i^), then g is called r~reducible modulo F. If \'ai{g) has an divisor in 
lm(F), polynomial g is called r~head reducible modulo F. If for any f £ F and for 
any monomial u occurring in / there is no r-divisors of u in lm(F) \ {lm{f )}, then F 
is called r-autoreduced. 

Definition 2.6. Similarly to the conventional division in Definition 2.2, an r—head 
normal form HNFr{h, F) and r— normal form NFr{h,F) of /i G M \ {0} modulo 
-F C K \ {0}, for a given restricted division r, is the polynomial h given by 



Here lm(/) |,. lm(/) Wy if ofy 7^ 0, lm(my/,) ^ lm(/i). And, respectively, \mh and 
all monomials in h are r— irreducible modulo lm(F). 

Definition 2.7. Given an ideal X c M, an order y- and a restricted division r, a finite 
subset G C M is called r— basis of T if 



As well as in the case of the conventional monomial division, r— basis G can also be 
defined by the relation 



(V/e/)[7VF(/,G) = 0]. 



(2) 





(V/eX)(3geG) [lm(g) |. lm(/) ]. 



(V/e/) [iVF,(/,G) =0]. 



From Definition 2.4 it follows that a r— basis is always a Grobner basis. 
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3. Involutive Division and Involutive Bases - Definition 

A natural way to introduce a restricted monomial division r is to indicate a certain subset 
X{u) C X of variables for a monomial m G M and to define for w G M 

where Mx(u) is the monoid of power products constructed from the variables in X{u). 
Definitions 2.5 and 2.6 deal with r-divisors taken from a fixed finite monomial set. By 
this reason it suffices to define an r-division for an arbitrary finite set of divisors. The be- 
low definitions, taken from 1 1 1 1, allowed to introduce a wide class of restricted divisions 
providing an algorithmic way for construction of r— bases 1 14 17 18J. 

Definition 3.1. We say that an involutive division L is defined on M if for any nonempty 
finite monomial set J7 C M and for any u eU there is defined a subset Mc{u, [/) C X 
of variables generating monoid C{u, U) = ^Mc{u,u) such that the following conditions 
hold 

1. u e C/ A uC{u, U) n vC{v, U)i=$^ue vC{v, U) \/ V e uC(u, U). 

2. veU A ve uC{u, U) => C{v, U) C C{u, U). 

3. uG^A V QU =^ C{u, U) C C{u, V). 

Variables in Mc{u, U) are called {£—) multiplicative for u and those in NMc{u, U) = 
X \ Mc{u, U) are called {C—) non-multiplicative for u, respectively. \f w £ uC{u, U), 
then u is called C— (involutive) divisor of w. In this case we shall often write u \c w. 

The restricted division given by Definition 3.1 called involutive since it generalizes 
the main properties of three different separations of variables used by Janet |38|, 
Thomas |39| and Pommaret |40| in the algebraic analysis of partial differential equation 
systems based on their completion to involution. Those separations satisfy 1 11 1 proper- 
ties 1-3 in Definition 3.1 and, hence, generate involutive divisions. 

Now, as particular examples of involutive divisions, we present Janet and Pommaret 
divisions 1 1 1 1 which are used most. In addition, we consider the example of involutive 
division introduced in 1 17 1. In doing so, we indicate multiplicative or no n- multiplicative 
variables only. The remaining variables are to be considered as non-multiplicative or 
multiplicative, respectively. 

Example 3.1. ( Janet division ). For each 1 < i < n divide J7 £ M into groups labeled 
by non-negative integers di, . . . ,di 

[di, . . . ,di] = {u eU\dj = degj(u), 1 <j <i}. 

x\ is multipUcative for w G C/ if degj^ [u) — maxjdeg^ (v) \v e [/}. For i > 1 Xi is mul- 
tiplicative for u e [di, . . . , di-i] when degj(M) = max{degj(?;) | v e [di, . . . , c?i-i]}. 

Example 3.2. ( Pommaret division ). For with dk > variables 

Xj,j > fc are multiplicative. For w = 1 all the variables are multiplicative. 

Example 3.3. ( Lexicographically induced division ). A variable Xi is non-multiplicative 
for w e [/ if there is w e [/ such that v <Lex u and degi{u) < degi{v), where ^Lex 
denotes the pure lexicographical monomial order. 
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In the last example one can replace the lexicographical order with any other (admissible) 
order to obtain another involutive division 1171 . The separation of variables into multi- 
plicative (Af) and non-multiplicative (NM) for divisions in Examples 3.1-3.3 and for 
the set of three monomials in three variables shown in Table 1 . 



Table 1. Separations of variables for U = {x^xs, xiX2, xix^} 



Element 
in U 


Involutive division 


Janet 


Pommaret 


Lex. induced 


M 


NM 


M 


NM 


M 


NM 


xlxi 


Xi,X2,X3 




X3 


Xl , X2 


Xl 


X2,X3 


XlX2 




X\ 




Xl 


Xl,X2 


X3 


X x^ 


^3 


Xl,X2 


X2,X3 


Xl 


Xi,X2,X3 





As it was shown for the first time by Zharkov and Blinkov 1121 . completion of multi- 
variate polynomial systems to involution gives an algorithm for construction of Grobner 
bases. Zharkov and Blinkov used Pommaret separation of variables and thus constructed 
what is called now Pommaret bases. They proved termination of their completion al- 
gorithm for zero-dimensional ideals whereas for positive-dimensional ideals Pommaret 
bases may not exist (as finite sets), and, thus, the algorithm may not terminate. 

Termination properties of completion algorithms are determined by the underlying 
involutive divisions and can be studied by completion of monomial sets 1 11 1. The fol- 
lowing two definitions elucidate this important aspect of involutive algorithms. 

Definition 3.2. Set C{U) ~ U„gt/ uM is called the cone generated by finite monomial 
set J7 C M and set Cc{U) = U„g(7 u£{u, U) is called C— (involutive) cone of U . 

Definition 3.3. A finite monomial setU ^ U is called C— completion of set U if C{U) = 
Cc{U).lf 

Cc{U)=C{U) (3) 

set U is called C-complete or C-involutive. If every finite set U admits £— completion, 
then involutive division C is called Noetherian. 

Janet division and Lexicographically induced division are Noetherian fl 1"17V As to 
Pommaret division, it is non-Noetherian | 111 what can be also seen by explicit com- 
pletion of the monomial set in Table 1 . The completion can be performed either by the 
special monomial algorithms fl 1"161 or by the general polynomial algorithms described 
below. The resulting sets are given by: 

Janet: {x\xy,^xiX2TXix\Tx\x2\, 

Pommaret : {x\xztXiX2,xix\,x\x2, ■ ■ ■ ,x\^'^X2^ ■ ■ ■ ,x{^'^xs,, . . G N>o, 

Lex. induced : {x\xy,, a;ia;2, a;ia::3, 2:1X2X3}. 

Apart from Noetherianity, algorithmically "good" involutive divisions must be continu- 
ous and constructive in accordance with two definitions that foUow. 
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Definition 3.4. A set U is called locally C— involutive or locally C— complete if 



An involutive division C is called continuous if conditions (0} imply equality for any 
monomial set [/ G M. 

Unlike the condition (|3} of involutivity /completeness, the conditions @ of local involu- 
tivity/completeness admit an algorithmic verification and are fundamental for involutive 
algorithms. In paper [1 IJ we found a criterion which allows to check continuity of an 
involutive division. 

Given a monomial order >- and an involutive division £, the separation of variables 
for elements in a finite polynomial set C M \ {0} is defined in terms of the leading 
monomial set lm(F) . Now we can introduce the concept of involutive basis as an r— basis 
in Definition 2.6 specified for an involutive division. But for all that we follow the defi- 
nition in 1.1 IJ and demand for involutive bases are to be autoreduced in accordance with 
Definition 2.7. 

Definition 3.5. Let X C M be a nonzero ideal, C be an involutive division and >- be 
a monomial order Then a finite >C-autoreduced subset G C K. is called (C—)involutive 
basis of X if 



If division C is continuous, then conditions ^ are equivalent II II to the following invo- 
lutivity conditions 



( V/ e G ) ( Vx, e NMc{\mif),\m{G)) ) [ NFc[x, • /, G) = ] . (6) 



Here NFc{f, F) denotes the £-normal form of / modulo F according to Definition 2.6 
where r— division replaced by ^-division. 

Definition 3.6. The product of a polynomial by its non-multiplicative (multiplicative) 
variable is called non-multiplicative (multiplicative) prolongation of the polynomial. 

Involutivity conditions (|5Jl give an algorithmic characterization of involutive bases much 
like the algorithmic characterization of Grobner bases by S"— polynomials established 
in 1 1 1. One can say that involutive bases are characterized by £— reducibility to zero of 
all the non-multiplicative prolongations of elements in the basis. 

The basic idea behind an algorithmic construction of involutive bases is to check 
the involutivity conditions. If they are not satisfied, then the nonzero /^—normal forms 
are added in a certain order to the system until all the involutivity conditions satis- 
fied. With all this going on, an algorithmically "good" involutive division, in addition 
to its Noetherianity and continuity, should satisfy the constructivity properties II II in 
accordance to the definition: 

Definition 3.7. An involutive division C is called constructive if for any U cM.,u ^ U, 
Xi e NMciu, U) such that u ■ Xi ^ Cc{U) and 



(Vw e U) (Vx e NMc{u, U)) [ u • x e Cc{U) ] . 



(4) 



(V/eX)(3.geG)[lm(.g) U lm(/) ] . 



(5) 



g 
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(Vu e U) {Mxj G NMc{v, U)) {v-xj nu-x^) [v ■ xj e Cc[U) ] 

the condition (Vw £ Cc{U)) [u ■ Xi ^ wC{w, U U {w}) ] holds. 

Constructivity ensures that in the course of the algorithm there no needs to enlarge the 
intermediate basis with multiplicative prolongations and only non-multiplicative prolon- 
gations must be examined for the enlargement (completion). Note that all three divisions 
of Examples 3.1-3.3 are continuous and constructive II 11171 . 

As any i — basis, an involutive basis is a Grobner basis. £-reducibility implies the 
conventional reducibility (i.e. reducibility with respect to the conventional division). But 
the converse is not true in general. By this reason, an involutive basis is generally redun- 
dant as the Grobner one. Moreover, a monic reduced Grobner basis is uniquely defined 
by an ideal and a monomial order L1.13il whereas this is not true for involutive bases as 
shows the following simple bivariate example 111 II . 

Example 3.4. Consider ideal in Q[x, y] generated by F = {x^y — 1, xy^ — 1}. For the 
lexicographical order with x y y the polynomial sets 

Gi = {xy^ -^,xy- y'^,x -y,y^ - I}, 

G2 = {x^y -l,x^ - y^,xy^ - y, xy"^ -l,xy- y'^, x-y,y^ -1} 
are both Janet bases of Id(F). 

However, for a constructive division one can define fl41 a minimal involutive basis 
which, similarly to a reduced Grobner basis, being monic is uniquely defined by an ideal 
and a monomial order. If G is such a monic minimal involutive basis, then for any other 
monic involutive basis Gi the inclusion G C Gi holds. For Example 3.4. the minimal 
Janet basis is {x — y,y^ — 1} and coincides with the reduced Grobner basis. 

4. Involutive Bases — Construction 

In the rest of the paper we assume that the input involutive division in the below al- 
gorithms is Noetherian, continuous and constructive. Our goal is to construct minimal 
bases. Having this in mind, we shall often omit the word "minimal". 

First, we present the simplest version of an algorithm for constructing involutive 
polynomial bases and illustrate its work by the bivariate polynomial set from Example 
3.4. Then we describe an improved version of the algorithm. 

In the below algorithm InvolutiveBasis I the whole polynomial data are partitioned 
into two subsets G and Q. Set G contains a part of the intermediate basis. Another part of 
the intermediate basis is contained in set Q and includes also all the non-multiplicative 
prolongations of polynomials in G which must be examined in accordance to the involu- 
tivity conditions (|5Jl. 

At the initialization step of lines 1-3, an element from the input polynomial set F is 
chosen such that its leading monomial has no proper divisors (in the conventional, non- 
involutive sense) among the remaining elements in Im(i^). Similar choice of an element 
from Q is done in line 7. Before insertion of a new nonzero element h into G done in line 
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15, all elements g E G such that 1111(5) is ^ proper multiple of \m{h) are moved from G 
to Q in line 13. 

Such a choice in lines 1 and 7 together with the displacement step in line 13 provide 
correctness of the algorithm. If it terminates, then the output basis obviously satisfies 
the conditions (|6j since Q becomes the empty set by the while-condition in line 18. 
Minimality of the output basis can be proved |41 1 by exactly the same arguments as used 
in the proof of minimality for the algorithm in paper 1 14J . 

As to termination of algorithm InvolutiveBasis I, it follows from Noetherianity of 
division C. Indeed, there can be only finitely many cases II II when a nonzero element h 
obtained in line 9 and the polynomial p selected in line 7 is £— head reducible. Between 
such events, the set Im G of leading monomials in G can only be completed in the main 
loop 4- 1 8 by finitely many monomials which are not non-multiplicative prolongations of 
monomials in Im G. All other completions are performed just by these prolongations. By 
the Noetherianity of C, the last completion must terminate in finitely many steps II II . 



Algorithm: InvolutiveBasis 1(,F,^, C) 



Input: £ R \ {0}, a finite set; -<, an order; C, an involutive division 


Output: G, a minimal involutive basis of Id(i^) 


1 


choose f G F without proper divisors of lm(/) in lm(F) \ {lm(/)} 


2 


G ■■={/} 


3 


Q:^F\G 


4 


do 


5 


h:=0 


6 


while Q ^ and /i = do 


7 


choose p G Q without proper divisors of lm(p) in \m{Q) \ {lm(p)} 


8 


Q:=Q\ {p} 


9 


h := NormalForm(p, G,-<,C) 


10 


od 


11 


ith^O then 


12 


for all {geG \ lm(/i) C lm(g)} do 


13 


Q:=QU{g}; G:^G\{g} 


14 


od 


15 


G:^G(J{h} 


16 


Q:^QU{g-x\geG,xe NMc{lm{g), lm(G)) } 


17 


fl 


18 


od while Q ^ d) 


19 


return G 



It should be noted that the above algorithm is distinguished from the algorithm in 1161 
by the conditions of the choice and of the displacement used in lines 1, 3 and 7, respec- 
tively. In paper [161 the polynomials with the smallest leading monomials with respect 
to the order >- were chosen, and those with higher leading monomials then lm(/i) were 
displaced. The choice made in 1161 is apparently more restrictive than that in the above 
algorithm InvolutiveBasis I, and with more replacements than those done in line 1 3 of 
the algorithm. 
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Algorithm: NormalForm(p, G,<,C) 


Input: p eR \ {0}, a polynomial; G C M \ {0}, a finite set; 




^, an order; C, an involutive division 


Output: h — NFc{p, G), the £— normal form of p modulo G 


1 


h :— p 


2 


while h ^ Q and /i has a term t /^—reducible modulo lm(G) do 


3 


take g E G such that lm(g) \c t 


4 


h-.^h-g-t/ lt(.g) 


5 


od 


6 


return h 



Subalgorithm NormalForm invoked in line 9 of the InvolutiveBasis I algorithm 
computes £— normal form in the full accordance with Definition 2.6 specified for 
division. Its termination immediately follows from the fact that /^—reductions form 
a subset of the conventional reductions, and the last reduction sequence is always fi- 
nite cm. 

£— reducibility of polynomial h is checked in line 2 of algorithm NormalForm. 
The check consists in search for an £— divisor of \m{h) among elements in lm(G). 
Polynomial set G, as constructed in the course of the main algorithm InvolutiveBasis 
I, is £— autoreduced at every step of the completion procedure I14I41I . Therefore, by 
property 2 in Definition 3.1, \m{h) may have at most one £— divisor. Hence, polynomial 
g satisfying the condition taken in line 3 is unique in G as an £— reductor of the term t 
in h. 

Now we illustrate algorithm InvolutiveBasis I by Example 3.4. In Table 2 we show 
the intermediate values of sets G in the 2nd column, Janet non-multiplicative variables 
NMj for elements in G in the 3rd column and elements in Q in the 4th column. Rows of 
the table contain these values obtained at the initialization and after every iteration of the 
main loop 4-18. In this case in lines 3 and 7 we selected the lexicographically smallest 
elements. 



Table 2. Computation of Janet basis for Example 3.4 



Steps of 


Sets G and Q 


algorithm 


elements in G 


NAIj 


Q 


initialization 


xy^ — 1 




Wy-1} 


iteration 


x^y — 1 








xy^ — 1 


X 


{x^y'^ — x} 




x-y 




{xy^ — lyx'^y — 1} 




x-y 








y3 -1 


X 


{x-'y- l,xy'-' -x} 




x-y 










X 


{} 
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The grave practical disadvantage of the presented algorithm is that it treats useless re- 
peated prolongations and does not use any criteria to avoid unnecessary reductions. Be- 
low we describe an improved version of involutive algorithm where unnecessary repeated 
prolongations are avoided and where the involutive analogues of Buchberger's criteria 
are enabled. For these purposes we need the next definition. 

Definition 4.1. An ancestor of a polynomial / e F C M \ {0} is a polynomial g E F 
of the smallest deg{lm{g)) among those satisfying f = g ■ u modulo ld{F \ {/}) with 
M € M. If deg(lm(5)) < dcg(lm(/)) (u ^ 1) the ancestor 5 of / is called proper. If an 
intermediate polynomial h that arose in the course of a completion algorithm has a proper 
ancestor g, then h has been obtained from g via a sequence of £-head irreducible non- 
multiplicative prolongations. For the ancestor g itself the equality \m{ai\c{g)) — lm(g) 
holds. 

Let now every element / G F in the intermediate set of polynomials be endowed 
(cf. (16|) with the triple structure 

P = {/, 5, vars} 

where 

pol(p) = / is the polynomial / itself, 
anc(p) =g is a polynomial ancestor of / in F, 
nmp(p) = vars is a (possibly empty) subset of variables. 

The set vars associated with polynomial / accumulates those non-multiplicative vari- 
ables of / have been already used in the algorithm for construction of non-multiplicative 
prolongations. It keeps information on non-multiplicative prolongations of polynomial / 
that have been already examined in the course of completion and serves to avoid useless 
repeated prolongations. 

Knowledge of an ancestor of / in F helps to avoid some unnecessary reductions by 
applying the involutive analogues of Buchberger's criteria described below. An improved 
version of algorithm InvolutiveBasis I named InvolutiveBasis II is given as follows. 
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Algorithm: InvolutiveBasis II(F. -<.L) 



Input: F € R \ {0}, a finite set; -<, an order; £, an involutive division 


Output: G, a minimal involutive basis of Id(F) or a reduced Grobner basis 


1 


choose / e -F without proper divisors of lm(/) in lm(F) \ {lm(/)} 


2 


T:={/,/,0} 


3 


Q:={{9,g,0}UGi^\{/}} 


4 


Q :=HeadReduce((3, T, -<, £) 


5 


while g 7^ do 


6 


choose p G Q without proper divisors of lm(pol(/)) in 




lm(pol(Q)) \ {lm(pol(/))} 


7 


Q.= Q\ M 


8 


if pol(p) = anc(p) then 


9 


for all {qeT\ lm(pol(g)) □ lm(pol(p))} do 


10 


Q:=QU{q}; T:=T\{q} 


11 


od 


12 


fi 


13 


h :=TailNormalForm(j3, T, ^, £) 


14 


T := T U {/i, anc(p), nmp(p)} 


15 


for all g G T and x e NMc{q, T) \ nmp(g) do 


16 


g := Q U {{pol(5) • X, anc(g), 0}} 


17 


nmp(g) := nmp(g) n NMc{q,T) U {a;} 


18 


od 


19 


g :=HeadReduce(g, T, ^, £) 


20 


od 


21 


return {pol(/) | / G T} or {pol(/) | / e T | / = anc(/)} 



Here and in the below algorithms, where no confusion can arise, we simply refer to the 
triple set T as the second argument in NMc, NFc, and HNFc instead of the poly- 
nomial set {g ~ pol(t) I f G T }. Sometimes we also refer to the triple p instead of 
pol(p). Besides, when we speak on reduction of triple set Q modulo triple set T we mean 
reduction of the polynomial set { / = \)o\{q) | g G g } modulo {g = pol(f) | f G T }. 

Apart from providing subsets Q and T of the intermediate basis with the triple struc- 
ture, the improved version contains extra lines 4 and 19. Here C— head reduction is done 
for the basis elements in Q modulo those in T. Then the remaining tail reduction is per- 
formed in line 1 3 to obtain the (full) £-involutive normal form. Furthermore, unlike the 
previous algorithm and due to presence of the third elements in triples, the set Q is en- 
larged in line 16 only with those non-multipUcative prolongations which have not been 
examined yet. 

In doing so, a new prolongation of a polynomial in T is inserted into Q with the 
ancestor of the polynomial. In the next Une the selected non-multiplicative variable x is 
added to the set of non-multiplicative variables already used. The intersection in line 17 
takes into account that some of these variables may turn into multiplicative through the 
contraction of T in line 10 and by virtue of the relation 3 in Definition 3.1. 
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The subalgorithm HeadReduce invoked in lines 4 and 19 of the main algorithm 
returns set Q which, if nonempty, contains the part of intermediate basis £— head reduced 
modulo T. 



Algorithm: HeadReduce (Q, T,^,£) 



Input: Q and T, sets of triples; ^, an order; C, an involutive division 


Output: £— head reduced set Q modulo T 


1 


a 


2 




3 


wniie o y= w uo 


4 


cnoose p o 


5 


S:^S\{p} 


6 


h ^HeadNormalForm T C\ 


7 


ith^O then 


8 


if lm(pol(p)) ^ lm(/i) then 


9 


Q Q U {/i, /i, 0} 


10 


else 


11 


g := Q U M 


12 


fi 


13 


else 


14 


if lm(pol(p)) — lm(anc(p)) then 


15 


for all {q e S \ anc(g) = pol(p)} do 


16 


S:=S\{q} 


17 


od 


18 


fi 


19 


fi 


20 


od 


21 


return Q 



Its own subalgorithm HeadNormalForm invoked in line 6 just computes the head 
normal form h of the polynomial pol(p) in the input set Q assigned to the set S at the 
initialization step (lines 1-2). If ft, ^ when lm(pol(p)) is £— reducible what is verified 
in line 8, then lm(ft,) does not belong to the initial ideal generated by { lm(pol(/)) | / G 
Q UT} 1141 . In this case the triple {h, h, 0} for h is inserted (Une 9) into the output set 
Q. Otherwise, the output set Q retains the triple p as it is in the input. 

In the case when h — 0, whereas pol(p) has no proper ancestors what is verified 
in line 14, all the descendant triples for p, if any, are deleted from S in line 16. Such 
descendants cannot occur in T owing to the choice conditions in lines 1, 6 and to the 
displacement condition of line 9 in the main algorithm InvolutiveBasis II. Step 14-18 
serves for the memory saving and can be ignored if the memory restrictions are not very 
critical for a given problem. In this case all those descendants will be casted away by the 
criteria checked in the below algorithm HeadNormalForm. 

The next algorithm HeadNormalForm, after initialization in lines 1-2, starts with 
verification in line 3 of £— head reducibility of the input polynomial h modulo polyno- 
mial set G := { pol(5) | 5 G T }. This verification consists in search for £— divisor of 



14 



V.P. Gerdt / Involutive Algorithms 



\m{h) in lm(G). If there is no such divisor the algorithm returns h in line 4. Otherwise, 
in the course of the search, the polynomial g E G is found such that lm(g) \c lm{h). 



Algorithm: HeadNormalForm(p, T,^,£) 



Input: T, a set of triples; p, a triple such that pol(p) = pol(/) • x, 




/ G r, a; G NMc{f, T); an order; £, an involutive division 


Output: h = HNFc{p, T), the £— head normal form of pol(p) modulo T 


1 


/i pol(p) 


9 
Z 




3 


if lm(/i) is £-irreducible modulo G then 


4 


return /i 


5 


else 


6 


take g eT such that lm(pol(5)) £ lm(/i) 


7 


if lm(/i) 7^ lm(anc(p)) then 


8 


if Criteria{p, g) then 


9 


return 


10 


fi 


11 


else 


12 


while ft- 7^ and lm(/i) is £— reducible modulo G do 


13 


take g in T such that \m{q) \c lm(/i) 


14 


h:^h-q-\t{h)/\t{q) 


15 


od 


16 


fi 


17 


fi return h 



If the set G is £— autoreduced, and this just takes place when algorithm HeadNormal- 
Form is invoked in line 6 of algorithm HeadReduce, then there is the only one such g 
in G. Similarly, the check whether the further £— reducibility takes place is done in line 
12, and if that is the case, then the corresponding unique reductor is taken in line 13. 

For the £-head reducible input polynomial pol(p) what is checked in line 3, the 
following four criteria are verified in line 8 

Criteria{p, g)=Gi{p, g) V C'2(p, g) V Gz{p, g) V Gi{p, g), 

where 

Gi{p,g) is true <==^ lm(anc(p)) • \m{a.nc{g)) — lm(pol(p)), 

G2{p,g) is true <==^ lcm(lm(anc(p)), lm(aiic(g))) C lm(pol(p)), 

C3 {p, g) is true 3 t <E T such that 

lcm(lm(pol(t)), lm(anc(p))) C lcm(lm(anc(p)), lm(anc((7))) A 
lcm(lm(pol(t)), lm(anc((7))) □ lcm(lm(anc(p)), lm(anc(5))), 

G4(p,.g)istrue4=» 3i eT Ay e N Mc{t,T) v/ithlm{pol{t))-y = lm(pol(p)), 
lcm(lm(anc(p)),lm(anc(t))) C lm(pol(p)) A idx(t,r) < idx(/, T), 

where idx{t, T) enumerates position of triple t in set T. 

Criterion Gi is Buchberger's co-prime criterion [13] in its involutive form P161. It is 
easy to see that under the condition lm(pol(5)) |£ lm(/i) of line 6 
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lm(anc(/)) • lm(anc(g)) — lm(pol(/)) lm(anc(/)) • lm(anc(5)) | lni(pol(/)), 

and in [16] the right-hand side form of the criterion was given. Criterion C2 was derived 
in lUl I as a consequence of Buchberger's chain criterion 1 13 1. Criteria C3 and C4 derived 
in paper 1 32 1 complement criterion C2 to the full equivalence to the chain criterion. Thus, 
the four criteria in the aggregate are equivalent 1 32 1 to Buchberger's criteria adapted to 
the involutive completion procedure. 

The last subalgorithm TailNormalForm completes the £— head reduction by per- 
forming the involutive tail reduction and is invoked in line 1 3 of the main algorithm In- 
volutiveBasis II. In that way, it returns the £— normal form of the input £— head reduced 
polynomial as given by Definition 2.6. 



Algorithm: TailNormalForm (p, T,-<,C) 



Input: p, a triple p such that pol(p) = HNFc{p, T); 




T, a set of triples; ^, an order; C, an involutive division 


Output: h ~ NFc{p, T), the £— normal form of pol(p) modulo T 


1 


G:={pol(5) IgeT} 


2 


h := pol(p) 


3 


while h has a term t £— reducible modulo G do 


4 


take g & G such that \m{g) \c t 


5 


h:^h-g-t/\i{g) 


6 


od 


7 


return h 



The main algorithm can output either involutive or reduced Grobner basis (or both) de- 
pending on the instruction used in line 21. It is easy to verify that the main algorithm 
together with its subalgorithms ensures that every element in the output involutive basis 
has one and only one ancestor This ancestor is apparently irreducible, in the Grobner 
sense, by other elements in the basis. Thereby, those elements in the involutive basis that 
have no proper ancestors constitute the reduced Grobner basis. 

Though involutive bases are usually redundant as Grobner ones, this specific redun- 
dancy, can provide more accessibility to information on polynomial ideals and mod- 
ules II21I33I . For example, unlike reduced Grobner bases, involutive bases give explicit 
simple formulae for both the Hilbert function |30| and Hilbert polynomial 1 15 16 1 of a 
polynomial ideal I. If G is an involutive basis of I, then the (affine) Hilbert function 
HFx{s) and the Hilbert polynomial HPx{s) are expressed in terms of binomial coeffi- 
cients as follows. 



^ ^ MGlm(G) ^ ^ 
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where is the number of multiplicative variables for u. The first term in the right 
hand sides of these formulae is the total number of monomials in M of degree < s. The 
sum in the expression for HFj{s) counts [301 the number of monomials of degree < s 
in the £— cone (see Definition 3.2) C£(lm(G)). This number coincides with the number 
of such monomials in the cone C(lm(G)) by the completeness condition C£(lm(G')) — 
C(lm(G)) for lm(G) that follows from Definition 3.5. 

There is a number of other useful applications of involutive bases in commutative 
and noncommutative algebra (see II20I21I33I and references therein). 

We conclude this section by observing that if one applies algorithm InvolutiveBasis 
II to construction Janet basis for the polynomial set in Example 3.4, then the interme- 
diate polynomial data coincide with those in Table 2, where G = { pol(p) | p £ T } 
denotes the T— part of the intermediate basis and Q stands for its Q— part. The only dis- 
tinction from functioning algorithm InvolutiveBasis I for this particular example is that 
prolongation xy'^ — x reduces to zero by criterion Gi. 



5. Efficiency Issues 

Computer experiments presented on the Web cite |http : / / invo . j inr . ru reveal su- 
periority (for most of benchmarks) of the CIC++ code 11161 implementing the Involutive- 
Basis II algorithm over the best present-day implementations of Buchnerger's algorithm. 

In this section we consider some efficiency aspects of algorithm InvolutiveBa- 
sis II as compared with the Buchberger's algorithm II13I when one treats all possible 
5— polynomials as critical pairs, under application of the both Buchberger criteria, and 
performs reductions by means of the conventional monomial division. 



5.1. Automatic avoidance of some useless critical pairs 



In the involutive case, whenever the leading monomial of the intermediate polynomial 
/, under consideration in the course of a completion algorithm, is £— reducible modulo 
another intermediate polynomial g, the corresponding S*— polynomial whose £— normal 
form is to be computed has the structure S{f, g) — f — gu. Monomial u G M con- 
tains only multiplicative variables for g. In terms of the ancestors / and g for / and g, 
respectively, this S*— polynomial corresponds to the conventional 5— polynomial 

lm(/) lm(5) 

If ancestor / is proper, then / has been obtained from / by a sequence of £— head 
irreducible non-multiplicative prolongations. Insertion of these prolongations into the 
intermediate basis can only be accompanied with tail reductions. 

Thus, in the involutive procedure only S*— polynomials (critical pais) of form Q 
must be processed by computing their involutive normal form. In this way some useless 
5— polynomials are automatically ignored in the course of involutive algorithms. 

As a simple example, consider the polynomial set presented in Table 3. The first 
column contains the polynomials. Their Janet non-multiplicative variables for a; ^ z 
and the corresponding prolongations shown in the second and third column, respectively. 
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Table 3. Avoidance of useless critical pair in involutive completion 



Polynomial 


NMj 


Prolongation 


5— polynomial 


pi=xy-l 








P2 = xz - 1 


y 




•S'(p2,pi) = yp2 - zpi 


PS =yz-l 


X 


XP3 


S{P3,Pi) =xp'i - zpi 



The forth column contains S*— polynomials to be examined in a completion procedure. 
We see that S*— polynomial S{p2,P3) is not considered at all. To avoid useless reduction 
of S{p2,P3), Buchberger's algorithm needs the chain criterion. 

5.2. Weakened role of criteria 

The fact that a number of useless S'— polynomials is automatically avoided in the in- 
volutive approach contributes in its much weaker dependence on the use of criteria as 
compared with Buchberger's algorithm. Already the first implementation of the involu- 
tive completion procedure for zero-dimensional ideals and Pommaret division done by 
Zharkov and Blinkov in Reduce 1 12 1 revealed higher computation efficiency of this pro- 
cedure in comparison with the Reduce implementation of Buchberger's algorithm. And 
at that time it was unexpected, since no criteria were used in II12I whereas Buchberger's 
criteria had been implemented in the Groebner package of Reduce. It is well-known that 
without the criteria Buchberger's algorithm becomes impractical even for rather small 
problems. 

It should be also said that our implementation in 1 16| of algorithm InvolutiveBasis 
II for Janet division, as presented on the Web site http : / / invo . j inr . ru exploited 
only criteria Ci and C2 (Sect.4). The effect of other two criteria C3 and C4 that are due 
to Apel and Hemmecke 1 32 1 (see also 1 35 1) for our algorithm is now under experimental 
investigation. 



Table 4. Influence of criteria 



Example 


Applicability 




Timin 


g (sec.) 




Ci 


C2 


C3 




Without 


Cl ^C2 




Ci Ci 


Cyclic6 


98 


2 


4 




0.14 


0.09 


0.11 


0.12 


Cyclic7 


698 


190 


22 




73.01 


46.54 


47.15 


58.72 


KatsuraS 


173 


1 


1 




29.54 


22.44 


22.25 


27.48 


Katsura9 


344 




1 




358.34 


282.25 


278.25 


337.52 


Colin3 




114 


169 


1 


736.15 


81.87 


55.10 


76.72 


Assur44 


89 


60 


171 


3 


10.76 


10.13 


8.83 


10.35 


Reimer6 


63 


235 


179 


12 


30.65 


17.06 


7.46 


9.69 



Table 4 contains data of computer experiments with the criteria for some benchmarks 
from 113 613 7 1 and for computation over the ring of integers. Columns 2-5 in Table 4 
show how many times the criteria are apphed in our algorithm for the degree-reverse- 
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lexicographical order. In so doing, criterion Ci is checked first, then sequentially C2, C3 
and C4 2. 

The last four columns of Table 4 illustrate the effect of the criteria for the C code [161 
running on an Opteron-242 computer under Gentoo Linux. Comparison of the timings 
with and without criteria shows that the use of criteria is not so critical for algorithm 
InvolutiveBasis II as for Buchberger's algorithm. For the latter problems of the Cyclic? 
size and larger become intractable if no criteria applied. 

Note, that the check of criterion C4 is more expensive from the computational point 
of view [351. That is why for all the examples in the table thsi criterion is adversely 
affected on the timings as Table 4 indicates. 

In any event, however, criteria C3 and C4 comprise qualitative improvement of the 
involutive algorithm and their quantitative effect has to be further studied. 

5.3. Smooth growth of intermediate coefficients 

We performed experimental investigation of intermediate coefficient growth for polyno- 
mial systems with integer coefficients as they are completed to involution by algorithm 
InvolutiveBasis II for Janet division. Our observation is that this growth is much more 
smooth than in the case of Buchberger's algorithm. Swell of intermediate coefficients 
is a well-known difficulty of this algorithm. Even in the case when coefficients in the 
initial polynomial set and in the Grobner basis are small, intermediate coefficients can 
be huge and lead to a dramatic slowing down of computation or to running out of a com- 
puter memory. The following example taken from 1421 nicely illustrates this behavior of 
Buchberger's algorithm 

Example 5.1. Consider ideal I = Id{F) in Q[x, y, z] generated by the polynomial set: 



Its Grobner basis for the degree-reverse-lexicographical order with x y y y z is small 
G — {x, 4 y'^+1, z^} whereas in the course of Buchberger's algorithm, as it implemented 
in Macaulay 2, there arise intermediate coefficients with about 80,000 digits |42|. As to 
algorithm InvolutiveBasis II, it outputs Grobner basis G or Janet basis GU{y z^,y^z^}, 
depending on the instruction in line 2 1 , with not more than 400 digits in the intermediate 
coefficients. For comparison with the timings for this example given in |42|, we run our 
code 1 16 1 on a 500 Mhz PC with 256 Mb RAM and computed the Grobner basis in 0.1 
seconds. 



For most of benchmarks from f3o37i and the degree-reverse-lexicographic order the 
intermediate coefficients arising in the course of the involutive algorithm grow up only in 
several times in comparison with coefficients in the output basis. For the same collection 
of examples as in Table. 4, the maximal lengths of the input, intermediate and output 

^The numbers in columns 2-5 of Table 4 for Cyclic6 widely differ from those in 1351 where another involutive 
algorithm was experimentally analyzed. 



F = 



8 x'^y'^ + 5xy^ + 3x^z + x^y z, 
x^ + 2 y^z'^ + I3y^ z^ + by z'^, 
8x^ + 12y^ + xz^ + 3, 
Ix'^y* + ISxy^z'^ + y^z^. 



(8) 
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coefficients, measured by number of 64 bit words occupied, are accumulated in columns 
2-4 of Table 5. Here we applied the first three criteria without the forth one. The last 
column shows the ratio of the length in the third and forth columns called "swell factor". 



Table 5. Coefficient size in 64 bit words 



Example 


Input 


Intermediate 


Output 


Swell factor 


Cyclic6 


1 


3 


1 


3.00 


Cyclic? 


1 


11 


5 


2.20 


KatsuraS 


1 


5 


4 


1.25 


Katsura9 


1 


8 


6 


1.33 


Cohn3 


1 


168 


19 


8.84 


Assur44 


1 


93 


19 


4.89 


Reimer6 


1 


4 


4 


1.00 



Only for a few examples among those we have already analyzed, the swell factor reached 
several dozens. For instance, for the "f855" benchmark it is 30. 

To our opinion, there are several peculiarities of the involutive completion procedure 
that may provide such a smooth behavior of the intermediate arithmetics. 

• Both selection of prolongations from Q (and the critical pairs of form among 
them) for the reduction process and the reduction process itself are more restric- 
tive in the involutive approach than in Buchberger's algorithm. For the last algo- 
rithm, as a result of numerous computer experiments, some strategies such as nor- 
mal strategy fTj] and mainly the "sugar" f43| we found as heuristically "good" 
restrictions for selection of critical pairs. In the above involutive algorithms the 
leading term of a prolongation selected for the tail reduction and insertion into T 
must not divide the leading terms of other prolongations. However, this restric- 
tion is far yet from fixing the prolongation to be selected. For the degree-reverse- 
lexicographic order, in our implementation II16I we select that with the minimal 
total degree of the leading monomial. By analogy with strategy in [131, this strat- 
egy can be considered as normal. 

In doing so, the £— head pre-reduction of Q modulo T done before selection of 
a prolongation for insertion into T is also an ingredient of the selection proce- 
dure. It should be noted that in the HeadReduce algorithm one can also do the 
£-head reduction only for a part of elements in Q. Then selection in line 6 of 
the main algorithm must be done from the head reduced part. However, in the 
case of such a partial £-head reduction of Q, the polynomial pol(p) m p ^ S 
chosen in line 4 of algorithm HeadReduce must not have proper divisors in 
lm({ pol((7) I g e }) \ {pol(p)}. The option to perform partial head reduction 
gives rise to a certain freedom in the selection strategy. In our implementation 1 161 
for the degree-reverse-lexicographic order we select for the head reduction all the 
elements in Q containing polynomials of the minimal total degree. 
As regards the reduction process itself, it is entirely fixed in the above involutive 
algorithms, since there can be at most one elementary involutive reduction for any 
term. This is because of uniqueness of an involutive divisor among the leading 
monomials of polynomial reductors in T (see Sect.4). 
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• The T-part of intermediate bases mostly contains some extra polynomials, in 
comparison with intermediate bases in Buchberger's algorithm that are redundant 
in the Grobner sense. These extra polynomials which are non-multiplicative pro- 
longations of their ancestors in T, actively participate in the reduction process and 
may also prevent intermediate coefficient swell as already noticed in 1301 . 

5.4. Fast search for involutive divisor 

As well as in any algorithm for constructing Grobner basis, the bulk of computing time 
for examples large enough is expended for reductions. In the involutive reduction process 
the most frequent operation is search for an involutive divisor for a given term among the 
leading terms of the intermediate polynomial set. For the InvolutiveBasis II algorithm 
this operation is performed in lines 3 and 12 of subalgorithm HeadNormalForm and 
in line 3 of subalgorithm TailNormalForm. These two subalgorithms, especially the 
former one, are invoked (for examples large enough) enormous number of times from the 
main algorithm and its subalgorithm HeadReduce. Thus, an optimal search for £-divisor 
is an essential ingredient of an efficient implementation of the involutive completion 
procedure. 

We have already stated that in the above algorithms at every elementary reduction 
step there is the only £-divisor Therefore, a wanted search is such that when there is a 
divisor it is located as quickly as possible. Otherwise, the search stops at the intermediate 
step as early as possible and signals that there is no divisor. 

For Janet division ^ we developed in 1 16| a special data structures called Janet trees 
which allow to organize the wanted search. A Janet tree is a binary search tree that 
takes proper account of properties of Janet division in Example 3.1 and whose leaves 
contain (pointers at) monomials among which the searching is done. Complexity bound 
for search in a Janet tree with the maximal total degree of the leaf monomials d in n 
variables \s 0{d + n). This is substantially lower than that for the binary search in a 
sorted set of monomials 1 16 1. 

We refer to L16J for more details on Janet trees, and illustrate their usefulness by an 
example. 

Example 5.2. Consider monomial set U — {x^y, xz, y^, yz, z^} and order x )~ y )~ z 
on the variables. According to this order which determines Janet separation of variables, 
assign number 1 to x, number 2 to y and number 3 to z. The below figure shows the 
structure of Janet tree for set U. 

The elements in U are located in the leaves. Each edge is associated with a certain vari- 
able, and each interior node has two integer indices. The root has 1 and as its indices. 
For other interior nodes the first index is the number of variable associated with the edge 
connecting the node with its parent. The second index is the degree of the variable such 
that, if the node is not the root, there exists a monomial in U which contains this variable 
in the degree indicated. In so doing, we consider the most compact form of the tree which 
is used by our C/C++ code. The first index of the right interior child is the number of the 
next variable. The left edge for the root is associated with the first variable. And the left 



Pommaret division admits sliglitly different trees. For extension to other divisions see im and m. 
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child of a node has always a higher second index (degree in the current variable) than its 
parent. 



Given a monomial u = x'^y^z^ {i, j, k € N> 0), we start in our search for a Janet divisor 
in U from the root. If z > we select x as the current variable with number 1, and 
compare i with the second index in the left child of the root ( 1 in our case ). 

If i = 1 we move to the child, and there is no further way to the left since the left 
child of our current node has degree (second index) 3 in the current variable x. Hence, 
the leaf monomials of the lest subtree cannot divide u. Thus, we have to look at the right 
child of our current node. If it exists, as in our case, we select the second variable y as 
the current one. Then we compare j with the second index (degree in y and in our 
case) of the right child. If j > we move to the right child. Since in our case the current 
node does not have left child, we change the current variable to z and compare k with 
the degree in z of the left child which is a leaf. If A: > 1 we go to the leaf and this means 
that the leaf monomial (xz in our case) is a Janet divisor we are looking for. Therefore 
for a monomial of the form xy'^z^ with j > 0, /c > 1 Janet divisor is xz. 

If « = 2 we cannot move to the left, since for the current variable x the second index 
of the left child is 3, and the leaf monomials in the left subtree (in our case the only 
monomial x?y) cannot contain Janet divisor of u. Therefore, a monomial of the form 
x^y^z^ has no Janet divisor in the tree. 

If i > 2 we go to the left child and then proceed with the next variable y, as the 
current one, by comparing j with the degree in y of the right child. It is leaf and has 
degree 1 . It follows that x^y is a Janet divisor for u with i>3,j>l,fc>0 whereas for 
i > 3, j = 0, A; > monomial u has no Janet divisor in U. 




(1,0) 



(1,3) 
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Similarly one can analyze the case i — when one has to go from the root to the 
right. 

This example illustrates the following general fact [161. For ^ given monomial, a 
Janet tree provides the unique path in the tree which either ends up with the leaf contain- 
ing the sought divisor or breaks in a certain interior node. The node is such that its left 
child, if any, has a higher degree in the current variable and there is no way to the right 
since the right child has a higher degree (in the next variable according to the order) than 
the given monomial. 

5.5. Parallelism 

Algorithm InvolutiveBasis II admits a natural and efficient parallelism. Naturalness is 
apparent from lines 4 and 19 of the main algorithm and from the structure of subalgo- 
rithm HeadReduce. The £— head reduction of polynomials in Q, that is, the most time- 
consuming part of the completion procedure can be done in parallel. As to efficiency, in 
our recent paper 1145 1 a slightly modified version of the algorithm oriented to the multi- 
thread computation was experimentally studied on a two processor Pentium III 700 Mhz 
computer running under Gentoo Linux. Some of experimental data extracted from those 
obtained in }A5\ and related to the benchmarks akeady used above, presented in Table 
6. The second and third columns show the timings for the one- and three-thread modes 



Table 6. Timings (in seconds) and speedup due to parallelism 



Example 


1 Thread 


3 Threads 


Speedup 


*1 th/H th 


Cyclic6 


0.79 


1.16 


-0.37 


0.68 


Cyclic? 


386.89 


182.86 


+294.03 


2.12 


KatsuraS 


119.92 


53.72 


+66.20 


2.23 


Katsura9 


1356.37 


587.82 


+768.55 


2.31 


Cohn3 


554.75 


222.69 


+332.06 


2.49 


Assur44 


73.93 


31.34 


+42.59 


2.36 


Reimer6 


88.99 


52.56 


+36.43 


1.69 



of parallel computation. The last two columns give absolute and relative speeding-up 
of the three-thread run in comparison with one-thread. Each of the threads was doing 
the involutive head reduction. We experimented with different numbers of treads, and it 
turned out that the maximal speed-up is achieved just for three threads. This is in confor- 
mity with the well-known observation that for multi-treading on a computer with SMP 
(Symmetric Multiprocessing) architecture an optimal number of threads is the number 
of processors exceeded by one. 

For such a small example as Cyclic6 the overheads of multi-threading surpass the 
computational effect of the parallelization, and we obtain some slowing down. For the 
other examples in Table 6 the speed-up is rather considerable. The relative speed-up 
greater than the number of processors can be explained by a better selection strategy dy- 
namically realizable in the multi-thread mode. In this case some head-reduced prolonga- 
tions may come into play earlier than in the one-thread mode and cause a faster chain of 
the successive reductions. 
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Since computational costs for Buchberger's algorithm are highly unstable with re- 
spect to the selection strategy for S"— polynomials, its experimental parallelization (see, 
for example, [461) does not reveal a reasonable gain from the parallelism. The involutive 
algorithm, according to our two-processor experiments, does not have this difficulty and 
rather stable to variation of the selection strategy. Moreover, if one performs in parallel 
the £— head reduction, then, as we observed in the three-thread run on the two-processor 
machine, it may even optimize selection strategy. At least, for a degree compatible term 
order, if the normal strategy is used, the parallelization apparently helps to select a pro- 
longation with the smallest total degree of the leading term. As we observed experimen- 
tally, in most cases this increases the speed of computation. 

Certainly, one has to run the algorithm on a machine with more number of processors 
to investigate its parallelization experimentally. Because of the result of two-processor 
benchmarking and the intrinsic parallel structure of the algorithm, one can expect its 
good experimental scalability. 

6. Conclusions 

As discussed above, experimental analysis of the described involutive algorithms imple- 
mented for Janet division shows that they form an efficient computational alternative to 
Buchberger's algorithm for construction of Grobner bases. 

Evidently, there are infinitely many different and algorithmically "good", i.e. 
Noetherian, continuous and constructive, involutive divisions |17| satisfying properties 
in Definition 3.1. However it is still an open question, if there is a better division in the- 
ory and/or heuristically better in practice than Janet division. In theory, one can say that 
division Ci is better than C2 if the inclusion 

(VC/ e M) (Vu e U) [ Mc, 2 Mc^ ] (9) 

holds for conventionally autoreduced monomial sets U, and there are sets U for which 
the inclusion is strict. In this case, by the same arguments as those used in |14| for 
the correctness proof, one can show that, for an identical selection strategy, algorithm 
InvolutiveBasis II with division Ci will process never more, but sometimes less number 
of non-multiplicative prolongations than with division £2- In addition, the corresponding 
output bases satisfy Gi C G2 with the strict inclusion for some ideals. 

In practice, running time depends not only on the number of prolongation treated and 
on the related size of the intermediate and final basis. Apart from other implementation 
aspects, running time depends also on costs of such frequently used and depending on the 
involutive division operations as re-computation of multiplicative and non-multiplicative 
variables under change of the intermediate polynomial set and search for involutive di- 
visors. However, it is fairly advisable to prefer (for implementation too) divisions which 
generate less number of non-multiplicative variables. 

To illustrate the last statement, confront Pommaret and Janet divisions for zero- 
dimensional ideals when Pommaret bases always exist. In such a case at first glance 
Pommaret division (Example 3.2) is favoured over Janet one (Example 3.1). Indeed, it 
does not require re-computing separation of variables for elements in the intermediate 
basis if other elements ai^e added or removed. Besides, Pommaret division admits trees 
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for searching involutive divisors that are very similar to Janet trees, and the search in 
a Pommaret tree is as fast as in the Janet tree with the same leaves. Nevertheless, as 
implemented in the form of algorithm InvolutiveBasis II, Pommaret division compares 
unfavourably with Janet division. We experimented with both divisions while creating 
the C code presented in L16J . and for all the benchmarks used Pommaret division lost in 
speed to Janet division. 

The reason is that Janet division surpasses Pommaret division in accordance to (|9j 
(see LI 1.47.1 ). Table 1 gives an example of the strict inclusion. To provide existence of 
the related Pommaret basis, U can be considered as the leading monomial set for zero- 
dimensional ideal, for instance, lA{x\x'i — l, xiX2 — 1, xiXg — 1). Since, as shown in 1471 . 
whenever Pommaret basis exists it coincides with (minimal) Janet basis, algorithm In- 
volutiveBasis II for Janet division produces the same output as for Pommaret division 
but with processing less number of prolongations. 

For the time being we do not know involutive divisions better in theory and/or heuris- 
tically in practice than Janet division. In the approach of Apel |30| involutive division is 
defined locally in terms of (admissible for) a monomial set as satisfying properties 1 and 
2 in Definition 3.1 and not necessarily property 3. In this approach, given a monomial set 
U, one can always construct a division which surpasses Janet division for U in the sense 
of relation (|9} and often the strict inclusion holds |30 33|. Based on this (local) concept 
of involutive division, Apel designed an algorithm for completion to involution with dy- 
namical construction of the best involutive division for every intermediate basis. This, 
theoretically very attractive, algorithmic procedure, is unlikely practical because of high 
computational costs of constructing the best intermediate divisions. A step forward to 
practicality of this approach was made by Hemmecke in the form of the SlicedDivision 
algorithm implemented in Aldor |34 35 1. Though this algorithm is far from being com- 
petitive in efficiency with algorithm InvolutiveBasis II, the both approaches to theory 
of involutive division are complementary |33| and their further development mutually 
helpful. 

One should also say that algorithm InvolutiveBasis II is apparently much less ef- 
ficient than the fastest modern-day algorithms designed by Faugere: F4 [48J recently 
built-in Magma (version 2-11) and Fc, |49|. These two algorithms are based on entirely 
different completion strategy and exploit the linear algebra methods for polynomial re- 
duction rather than elementary reduction chains used by Buchberger's algorithm and the 
involutive algorithm. We have good reason to think that the algorithmic ideas of Faugere 
can be incorporated into the involutive methods too. 
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